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ABSTRACT 

We investigate the metal enrichment of the intracluster medium (ICM) in the framework of 
hierarchical models of galaxy formation. We calculate the formation and evolution of galaxies 
and clusters using a semi-analytical model which includes the effects of flows of gas and 
metals both into and out of galaxies. For the first time in a semi-analytical model, we calculate 
the production of both a and iron-peak elements based on theoretical models for the lifetimes 
and ejecta of type la and type II supernovae (SNe la and SNe II). It is essential to include 
the long lifetimes of the SNIa progenitors in order to correctly model the evolution of the 
iron-peak elements. We find that if all stars form with an IMF similar to that found in the solar 
neighbourhood, then the metallicities of O, Mg, Si and Fe in the ICM are predicted to be 2-3 
times lower than observed values. In contrast, a model (also favoured on other grounds) in 
which stars formed in bursts triggered by galaxy mergers have a top-heavy IMF reproduces 
the observed ICM abundances of O, Mg, Si and Fe. The same model predicts ratios of ICM 
mass to total stellar luminosity in clusters which agree well with observations. According to 
our model, the bulk of the metals in clusters are produced by and brighter galaxies. We 
predict only mild evolution of [Fe/H] in the ICM with redshift out to z ^ 1, consistent with the 
sparse data available on high-z clusters. By contrast, the [O/Fe] ratio is predicted to gradually 
decrease with time because of the delayed production of iron compared with oxygen. We find 
that, at a given redshift, the scatter in global metallicity for clusters of a given mass is quite 
small, even though the formation histories of individual clusters show wide variations. The 
observed diversity in ICM metallicities may thus result from the range in metallicity gradients 
induced by the scatter in the assembly histories of clusters of galaxies. 

Key words: stars: luminosity function, mass function - galaxies: clusters: general - galaxies: 
formation - large-scale structure of the universe 



1 INTRODUCTION 

The bulk of the baryons in galaxy clusters are in their hot gas com- 
ponent, the intracluster medium (ICM), which typically contains 
5-10 times the mass found in cluster galaxies. The metal content 
of the ICM is sensitive to the star formation histories of the cluster 
galaxies and to the way in which metals are ejected from galaxies 
in supernova explosions and winds. In particular, the abundances 
of different species can be used to constrain the relative numbers of 
type la (SNe la) and type II (SNe II) supernovae that occurred in the 
cluster galaxies. This is because the abundances of the ct elements 
(e.g. O, Mg) are driven primarily by the rate of SNe II, whereas 
the production of iron is dominated by SNe la. The abundance of 
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a elements in the ICM is found, from X-ray measurements, to be 
comparable to abundances in the solar neighbourhood, whereas the 
abundance of iron is only about 30% of that found locally (e.g. 
.Mushotzkv et aL 199 6). This discrepancy has led several authors 
to propose that the initial mass function (IMF) of stars forming in 
cluster galaxies could be substantially different from that inferred 
for star formation in less dense environments (e.g. iRe nzini et alj 
'1993; Zeof & Silk 1996; Mushotzkv et al. 1996!; IWdarninllOOa 
iTornatore et aL.2004; Maoz & Gal-Yam„2004.) . 

The chemical evolution of galaxies has b een modelle d 
extensively in the monolithic collapse scenario (e.g.,|Larso3l96i). 
Simple infall models with delayed metal enrichment due to SNe la 
have been used to demonstrate that the metallicity distribution of 
stars in the solar neighbou rhood ca n be exp la ined by bottom-heavy 
IMFs as proposed by Salpetei] j l955t) j Matteucci & Gre^^id 
.1986; Matteucci & Francois! Il989l; iTsuiimoto et alJ Il99a 
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The metal enrichment of the ICM has also been studied in the 
framework of monolithic collapse, with the conclusion reached that 
Salpeter-like IMFs produce too few heavy elements to match the 
observed metal abundances in the ICM ( David, Forman & Jones 

[1991; Matteucci & Gibson 1995; Loewenstein & Mushotzkv 1996; 
Gibson & Matteuccii il997aih: .Moretti, Portinari & Chiosi .2003iV 
To reconcile the predictions with the observed ICM abundances, 
these studies concluded that a flatter IMF than Salpeter is required, 
such as the Arimoto-Yoshii IMF, with a; ~ 1 ( Ari moto & Yoshiii 
Il987h . where the number of stars per logarithmic interval in mass 
is given as dN /d\nm oc . (For comparison, the Salpeter IMF 
has X — 1. 35.) This conclusion has been confirme d by numerical 
simulations iValdam ini 2003; Torna tore et all2004t) . 

There is now overwhelming evidence supporting the gravi- 
tational insta bility paradigm for the formation of struc ture in the 
Universe (e.g. Peac ock et al. . l200lHBaugh et al. l200 4i^. Any real- 
istic model of galaxy formation should therefore be placed in the 
proper context of a model that follows the hierarchical formation 
of cosmic structures, su ch as the successful A-dominated cold dar k 
matter (ACDM) model jPercival et al. l2002tlSDergel et al. 120031) . 



Semi-analytic models of galaxy formation provide a framework 
within which the physics of galaxy formation can be followed at 
the same time as dark matter halos are built up through mergers and 



the accretion of smaller objects (e.s..'White & Frenklll991 


|Cole| 


199 it 


Kauffmann. White & Guiderdoni 19931; ICole et alj 


199i 


2000; Somerville c& Primack 1999; Nasashima et al."200r, 


2002; 


Menci et al. 2002; Hatton et al. 


2003). These models are able to 



reproduce a wide range of properties of the local galaxy popula- 
tion and to make predictions for the high redshift universe (e.g. 
[Saugh et al. 2004b). To date, the majority of semi-analytic models 
have used the instantaneous recycling approximation for chemical 
enrichment, and only considered metal production by SNe II. Such 
models have been used to investigate the stellar and gas metallici- 
ties of disk and elliptical galaxies and of damped Lya systems, with 
broad agreement obtained between model predictions and obser- 
vatio ns I Kauffmann 1996; Kauffmann & Chariot 1998; Cole et al. 
200Q; .Somerville P rimack & Fabet 2001; Nagashima & Yoshii 
2004 i lOkoshi et al 12 004). Semi-analytic models based on instanta- 
neous recycling have also been applied to enrichment of the ICM by 
Kauffm ann & Chariot (1998) and De Lucia. Kauffmann & White 
i2004h . as well as by Nagashima & Gouda (2001). The former two 
papers find that the observed ICM metallicities are only reproduced 
in models in which the chemical yield is twice the standard value 
for a solar neighbourhood IMF, in agreement with the conclusions 
from monolithic collapse models. 

The first attempt to take into account metal enrichme nt due to 
SNe la in hierarchical galaxy formation models was that of lThoma^ 
^99^ and Thomas & Kauffmann ( 1999). However, their calcula- 
tions used only the star formation histories for spiral and elliptical 
galaxies extracted from the Munich group's semi-analytical model, 
and assumed closed-box chemical evolution. They did not include 
the effects of gas inflows and outflows on the chemical enrichment 
histories of galaxies. The metal enrichment due to SNe la was con- 
sistently included in a s emi-analytical model for the first time by 
iNagashima & Okamotol feoQi), who took into account the recy- 
cling of metals between stars, the cold gas in the disk and the hot 
halo gas. While the lifetime of SNe la's was simply assumed to be 
a constant 1.5Gyr, their model provides good agreement with the 
observed distributions of [Fe/H] and of [O/Fe] vs. [Fe/H] for both 
solar neighbourhoo d and bul g e stars, based on a Salpeter-like IMF. 

Recently, Bau gh et al. ] <2004d) showed that the observed 
number counts of faint sub-mm sources can be matched in hier- 



archical models if starbursts triggered by galaxy mergers are as- 
sumed to take place with a top-heavy IMF (with slope x — 0). 
Since varying the slope of the IMF affects the abundance ratios of 
Q elements to iron, an important test of this model is to examine 
its predictions for the metallicity of the ICM. In this pape r, using 
essentially the same model as used bylBa ugh etZ1<2004lJ) . we m- 
vestigate the metal abundances of the ICM. We ad opt a standard, 
"bottom-heavy" IMF proposed by Kennicutt 1 1983"), with x = 1.5 
for m ^ IMq, for quiescent star formation in disks, and a top- 
heavy IMF with a; = in starbursts. While the IMF for quies- 
cent star formation is reasonably well determined from observa- 
tions in the solar neighbourhood (e.g. Krouoa 2002), the form of 
the IMF appropriate to starbursts is still unclear. There is, however, 
some tentative observational evidence suppo rting a top-heavy IMF 
jTsuru et alll997tlSmith & Gallaghej200ll) . 

In this paper, we incorporate metal enrichment due to SNe la 
as well as S Ne II into the GALFORM semi-analytic galaxy forma- 
tion model jCole et al . 2000), to see whether the model with a top- 
heavy IMF is able to reproduce the metal abundances of the ICM. 
The model includes both inflows of gas and metals into galaxies due 
to gas cooling in halos, and ejection of gas and metals by galactic 
winds. We calculate the production of different elements by SNe la 
and SNe II based on theoretical yields as a function of stellar mass, 
integrated over the model IMF. The lifetimes of SNe la are consis- 
tently calculated for the adopted stellar IMFs, while instantaneous 
recycling is assumed for SNe II. We compute ICM abundances for 
O, Mg, Si and Fe for clusters of different masses, and compare these 
with observational data. In making these comparisons, we correct 
the observed abundances for radial metallicity gradients; these cor- 
rections are important for Si and Fe. 

The paper is set out as follows. In Section 2, we briefly de- 
scribe the GALFORM model. In Section 3, we provide a detailed 
explanation of how we model metal enrichment. In Section 4, we 
show what our model predicts for the present-day luminosity func- 
tion of galaxies, for different assumptions about feedback and ther- 
mal conduction. In Section 5, we present our predictions for the 
abundances of different elements in the ICM, and for the ratios of 
ICM mass to stellar luminosity and to total mass of dark halos, 
as functions of cluster X-ray temperature, and compare with ob- 
servational data on present-day clusters. Section 6 compares the 
model with observational data on the evolution of ICM abundances 
with redshift. Section 7 examines the predicted metal enrichment 
histories of clusters in more detail. Finally, Section 8 presents our 
conclusions. The Appendix describes how we correct the observa- 
tional data for radial abundance gradients, to estimate global ICM 
metallicities. 

Throughout this paper, the cosmological parameters are fixed 
to be Ho = 0.3, r^A = 0.7, /i = Ho/100 km s"^Mpc"^ = 
0.7 and as = .93. The solar abundances used are as given by 
iGrevesse & Sauva l 1 1998). Som e observational data based on the 
old solar abundances given by lAnders & Grevessel 1 1989) have 
been corrected to be compatible with the new more recently de- 
termined solar values. 



2 THE BASIC GALAXY FORMATION MODEL 

We calculate galaxy formation in the framework of the C DM model 
using the semi-analytical model GALFORM described in lCole et alJ 
(2000, hereafter CLBF), with some modifications as described in 
Baugh et al. (2004b, hereafter B04). The semi-analytical model 
calculates various processes as follows: (a) dark matter halos form 
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by hierarchical merging; (b) gas falls into these halos and is shock- 
heated; (c) gas which cools in halos collapses to form rotationally 
supported galactic disks; (d) stars form quiescently in disks; (e) 
supernova explosions inject energy into the gas, ejecting some of 
it from the galaxy; (f) dynamical friction brings galaxies together 
within common dark halos, leading to galaxy mergers; (g) galaxy 
mergers can cause the transformation of stellar disks into spheroids, 
and can also trigger bursts of star formation. The model also com- 
putes the luminosities and colours of galaxies from their star for- 
mation histories and metallicities, using a stellar population model 
and including a treatment of dust extinction. 

The CLBF model assumed a cosmic baryon fraction 0,^ = 
0.02, that all stars formed with a Kennicutt (1983) IMF, and that 
gas ejected from g alaxies by supernova feedback was retained in 
the host dark halo. iBenson et alJ i2003ah showed that if the more 
recent estimate fib ~ 0.04 was used, the original model predicted 
too many very luminous galaxies. They proposed two possible so- 
lutions to this problem: either (i) thermal conduction suppresses the 
cooling of gas in high mass halos, or (ii) gas is ejected from the ha- 
los of galaxies by superwinds. B04 adopted the superwind solution, 
but showed that in order to explain the numbers of Lyman-break 
and sub-mm galaxies observed at high redshift, further changes are 
required in the model: they proposed changes in the star formation 
timescale in galaxy disks, in the triggering of starbursts by galaxy 
mergers, and in the IMF of stars formed in bursts. In this paper, we 
adopt the B04 version of our galaxy formation model as the stan- 
dard one, except that we retain both the thermal conduction and 
superwind variants. Our implementation of these two processes i s 
similar to, but somewhat simpler than that o f lBensonetai]<2003iJ) . 
For comparison purposes, we will also present some results for the 
original unmodified CLBF model. In the remainder of this section, 
we give some more details about some of the ingredients of our 
galaxy formation model. Our modelling of chemical enrichment is 
described in Section|3| 

Halo mergers and gas cooling: The merger histories of dark 
halos are generated for parent halos at z = 0, for a given power 
spectrum of density fluctuations, using a Monte Carlo scheme 
based on the extended Press-Schechter (PS) formalism (Press & 
Schechter 1974; CLBF). The lowest mass progenitor dark halos 
are assigned baryons in the form of diffuse gas only, according to 
the universal baryon fraction fib/fio- The diffuse gas is raised by 
shock heating to the virial temperature corresponding to the depth 
of the halo potential well. A fraction of the diffuse hot gas cools and 
falls into the centres of the host dark halos, where it accretes onto 
the disk of the central galaxy. The rate at which hot gas cools within 
a dark halo is computed by assuming a spherically symmetric gas 
distribution and using metallicity-dependent cooling functions. The 
predictions of this simple cooling model com pare favourably with 
the results of direct gas-dynamics simulations iYoshida et alJZOOa 
^Uy et al. 2003). 

Star formation in disks: The rate at which stars form from the 
cold gas in a galactic disk is given hy = A/coid / t, , where i\fcoid 
is the mass of cold gas. CLBF used a star formation timescale that 
depended on the dynamical time of the galaxy. Here, we adopt a 
prescription without any scaling with dynamical time, following 
B04: 



CLBF, in which cold gas is reheated and ejected into the halo by 
supernovae explosions according to the prescription: 



r° (V'disk/200 kms-i) 



(1) 



where t» and a* are free parameters independent of redshift, 
whose values are set by requiring that the model reproduce the ob- 
served gas mass to luminosity ratios of present-day disk galaxies. 
We employ the supernova feedback model described by 



where 

/3 = (Kiisk/Vhot)" 



(2) 



(3) 



and Viiot and Ohot are free parameters. The latter is fixed to Ohot = 
2 and Viiot is tabulated in TableQfor the models considered in this 
paper. 

As described above, for the standard baryon fraction em- 
ployed in the model, fit, = 0.04, our simple cooling and feedback 
models result in the formation of too many giant luminous galax- 
ies compared with recent determinations of the galaxy luminosity 
function (e.g. Norbere et al. 2002). We have investigated two dif- 
ferent approaches to solving this problem, which we now describe: 
(i) the suppression of gas cooling by thermal conduction and (ii) 
the suppression of star formation by galactic superwinds. 

Thermal conduction: In the first approach, thermal conduc- 
tion of energy from the outer to the inner parts of the gas halo acts 
to balance the energy removed from the gas by radiative cooling. 
This scenario is suggested by the lack of emission lines from low 
tempe rature gas in the X-ray spectra of clusters (e.g. lPeterson et all 
2001). To model this phenomenon in a simple way, we completely 
suppress cooling in halos whose circular velocity exceeds Kond- 
An estimate of the value of Vcond and its dependence on redshift 
can be derived by balancing the cooling rate with the heating rate 
due to thermal conduction: 



noniA(r) ^ V • (kVT) ~ k{T)T/W 



(4) 



where Uc and ni are the number densities of electrons and ions 
respectively, estimated here as their mean values within the virial 
radius, A(r) is the cooling function, k.{T) is the thermal conduc- 
tivity, Tj^_the_tem£erati^^ and R is the radius of the halo 
(e.g.. lTakahara & TakahM3ll98li) . Assuming a mean halo density 
given by the spherical collapse model, the Spitzer form of conduc- 
tivity without saturation, k cx T^/^, and a power law cooling rate 
for thermal bremsstrahlung, A ocT^''^we obtain 



Vcond = Kond(l + zf''^- 

Here V^ond is given by 



,1/2 

K"o„d ^ 200 ^ 
fsp 



Av 
100 



1/4 



km s 



(5) 



(6) 



where /hot is the ratio of hot gas to total baryon mass in 
a halo, fsp gives the ratio of the conductivity to the Spitzer 
value, and Av is the mean overdensity of a virialized halo rel- 
ative to the average density of the universe. (We have assumed 
A = 1.7 X 10~^^ergcm^s~^(r/K)^/^ and ksp = 1-0 x 
lO'^^erg cm^^s^^K"^ (T/K)^''^.) For a flat uni verse with presen t 
day density parameter fio = 0.3, Av ^ 300 jEke et alJll99^ . 
Assuming /hot ~ 0.8 (a realistic value in our models) then gives 
Kond — 240kms~^ for conduction at the Spitzer rate (/sp — 1). 
To match the bright end of the observed luminosity function, we 
show in Section 4 that we require 

Kond = lOOkms-^ this im- 
plies /sp — 20 if /hot = 0.8. While the cut-off circular velocity 
required increases to ~ 200 km s^^ if the redshift dependence of 
Vcond is neglected, the implied efficiency of conduction still seems 
unrealistically high. 

If conduction at this level is not viable, another mechanism 
needs to be invoked in order to prevent too much gas from cooling 
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in massive iialos. One plausible, but as yet relatively unexplored 
alternative mechanism is the heating of the hot gaseous halo b y 
emission from an active galactic nucleus re.g. lOranato et alJ2004) . 
This may be particularly relevant in the more modest mass ha- 
los in which the merger rate peaks jKauffmann & Haehnelt "2OO0I: 
lEnoki. Nagashima & Goud'3l2003h . Here we do not enter into the 
the details of the growth and fuelling of supermassive black holes 
in galaxy mergers. However, we retain the freedom to choose ap- 
parently implausibly low values of T4ond on the grounds that other 
phenomena in addition to conduction may be operating within the 
halo, with the same consequences, namely the suppression of gas 
cooling in halos above a certain mass. 

Superwinds: In the second approach we have taken to prevent 
the formation of too many bright galaxies, we assume that gas is 
ejected out of galaxy halos by superwinds. The superwind is driven 
by star formation, and ejects cold gas out of the galaxy and out of 
the halo at a rate given by 



where 

/3sw = 



: /sw min [1, (K/200 km s"^)"^] 



(7) 



(8) 



/sw is a parameter, and Vc is the circular velocity of the galac- 
tic disk for quiescent star formation or the circular velocity of the 
bulge in a starburst. In practise, we adopt the same values for the 
strength of the superwind, /sw, in the quiescent and burst modes 
of star f ormation. This is a sim plified version of the scheme pro- 
posed bv lBenson et alJ f2003ah . The form adopted for /sw means 
that the superwind becomes less effective at ejecting cold gas from 
galaxies with K > 200kms^^. The gas (and associated metals) ex- 
pelled from the halo are assumed to be recaptured once the halo is 
subsumed into a more massive object with a halo circular velocity 
in excess of Kocap- The free parameters in the superwind feedback 
are thus /sw and Kccap. 

Galaxy mergers and starbursts: Following a merger of dark 
matter halos, the most massive progenitor galaxy becomes the cen- 
tral galaxy in the new halo, and the remaining galaxies become 
its satellites. The orbits of the satellite galaxies decay due to dy- 
namical friction. If the dynamical friction timescale of a satellite 
is shorter than the lifetime of the new halo (defined as the time 
taken for its mass to double), the satellite will merge with the cen- 
tral galaxy. The consequences of the galaxy merger are determined 
by the mass ratio of the satellite to the primary, Msat/Mccn. If 
Afsat/Mccn > /ciiip, the merger is classified as a major merger; 
the merger remnant is a spheroid and any cold gas is consumed in a 
starburst. The duration of the starburst is proportional to the dynam- 
ical timescale of the bulge. For mergers with M^bx / Mcc-n < ,/ciiip, 
the merger is classified as a minor merger. In this case, the stars of 
the satellite are added to the bulge of the primary, but the stellar disk 
of the primary remains intact. If M^at / Mccn > /burst and the ratio 
of cold gas mass to total stellar and gas mass exceeds /gas.crit, the 
minor merger triggers a burst of star formation in which all the gas 
from the primary disk participates. We adopt the same values for 
the parameters governing the physics of galaxies mergers as used 
by B04: Uup = 0.3, /burst = 0.05 and /gas.crit = 0.75. We refer 
the reader to that paper for further discussion of the merger model. 

The standard model presented in this paper in both its ther- 
mal conduction and superwind variants is able to successfully re- 
produce many observed properties of present-day galaxies, such as 
luminosity functi ons, gas frac t ions, and sizes of galaxy disks, as 
will be shown in iLacev et'aP i2005l) . The comparison of the su- 
perwind variant of the model with observations of Lyman-break 



quiescent 
star formation 




recycling^ 
galaxy merger 



BULGE STARS 



Figure 1. Diagram showing the different baryonic components included in 
the model, and the mechanisms by which mass and metals are exchanged 
between them. 



and sub-mm galaxies at high redshift has already been published 
in Baugh et al. ( 2004b). The predictions for stell ar metallicities of 
galaxies will be presented in a subsequent paper iNagashima et alj 



3 STELLAR EVOLUTION AND METAL ENRICHMENT 

Our galaxy formation model thus has baryons in four different com- 
ponents: stars in galaxies, cold gas in galaxies, hot gas in halos, 
and hot gas outside halos. There are exchanges of material between 
stars and cold gas, cold gas and hot gas inside halos, cold gas and 
hot gas outside halos, and hot gas outside and inside halos, as de- 
scribed in the previous section. In addition, mergers of halos mix 
their reservoirs of hot gas, and mergers of galaxies mix the stel- 
lar and cold gas components of the two galaxies. The exchanges 
of mass and metals between the different baryonic components are 
shown schematically in Fig.Q Our model keeps track of the effects 
of all of these mass exchanges and merging events on the metal- 
licities of the different components. Finally, we must calculate the 
production of new heavy elements inside stars and their ejection 
back into the gas phase. We describe now how we model this last 
process. 



3.1 Evolution of a single stellar population 

Stars lose mass through stellar winds as they evolve. Some stars 
ultimately undergo supernova explosions which pollute the inter- 
stellar medium with newly produced metals. The rates of both 
of these phenomena depend upon the mass of the star. The to- 
tal amount of gas recycled and metals produced per unit stel- 
lar mass depends, therefore, on the adopted IMF. In this Sec- 
tion we explain how we obtain the gas restitution rate, the rates 
of SNe II and SNe la and the chemical yield for a single stel- 
lar population. This information is then incorporated directly into 
the GALFORM model, which predicts the star formation histories 
of individual galaxies, in order to obtain these quantities for the 
model galaxies. The discussion below is based upon that given by 
iLia. Portinari & Carraro' 1 2002) with slight differences. Stellar data 
are extracted from Portinari. Chiosi & Bressan ( 1998) for massive 
stars and Marigo (2001) for intermediat e and low mass stars. The 
chemical yields of SNe la are taken from lNomoto et ail il997h . 

We define the IMF as the number of stars per logarithmic in- 
terval of stellar mass per unit total mass of stars formed: 
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Figure 2. The initial mass functions (IMFs) of star formation used in tliis 
paper, 0(m) = dN/ d In m. For quiescent star formation we use the Ken- 
nicutt (1998) IMF (solid line) and for starbursts we use a flat IMF (dashed 
line). In both cases the IMFs are normalised to be unity when integrated 
over the whole range of mass. 



(m) = dA''/dlnm oc m 



(9) 



normalised to unity when integrated over the whole range of mass 



m 4'{m) d\nm 



(m) dm - 



(10) 



so that (/)( m) has units of M ^^. In this paper, we use the IMF pro- 
posed by iKennicutd il983h for quiescent star formation in disks, 
and a flat IMF in starbursts triggered by galaxy mergers. The 
[Kennicutt 1 1983) IMF has x = 0.4 for m < IMq and a; = 1.5 
for m ^ IMq. In starbursts, a; = 0. In both modes of star 
formation, the lower and upper stellar mass limits of the IMF 
are {md,mu) ~ (O.ISM©, I2OM0). Fig.|2|shows the quiescent 
(solid line) and starburst (dashed line) IMFs. 

In addition to the visible stars with m > O.ISM©, we allow 
for brown dwarfs with m < O.ISMq. The mass fraction of brown 
dwarfs is quantified by T, which is the ratio of the total mass in 
stars to the mass in visible stars for a newly formed population of 
stars. Thus, the mass fraction in visible stars is 1/T. CLBF adopted 
the Kennicutt IMF in both the quiescent and burst modes of star 
formation, and required T = 1.38 to match the predictions of their 
fiducial model to the observed present-day galaxy luminosity func- 
tion around L, . The other models we consider in this paper have 
T = 1. 

The recycled mass fraction from an evolving population of 
stars is a function of time. For a population of stars formed with an 
IMF (;/>(m) at time t — 0, the total mass fraction released by time t 
is 



[m - Afr(m)] (j){m) — , 



(11) 



M(t) 



where M{t) is the mass of a star reaching the end of its life at time 
t, and Mi{m) is the mass of the remnant left by a star with initial 
mass m. In our model, since the timestep we use is roughly com- 
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Figure 3. The cumulative returned mass fraction, t), from evolved 

stars after time t has elapsed from the epoch of star formation. The sohd 
and dashed lines indicate the E{i^ t) we use for the Kennicutt and top- 
heavy IMFs, respectively. Because we assume the instantaneous recycling 
for SNe II, the values are constant until r(8MQ), where r(m,) is the life- 
time of a star of mass m. The dotted lines show how these results change 
when the instantaneous recycling approximation is dropped. 



parable to the lifetime 4 x lO^yr of an %Mq star, the lowest mass 
star that results in a SNe II, we can justifiably use the instantaneous 
recycling approximation for SNe II. We therefore instantaneously 
release the mass recycled by stars with M 9>Mq. In Fig.|3| we 
show t) for both IMFs considered. The dotted lines show 

how these predictions change if the assumption of instantaneous 
recycling is relaxed. 

The cumulative number of SNe II explosions up to time t (per 
Mq of stars formed), -Rii(^ t), and the mass fraction of metals of 
the i-th element, Ei[^ t), released up to time t are given respec- 
tively by 



(j>{m) , 

max{M(t),8MQ} ^ 

' , , / N , / N dm 
Mi(m)(l){m) , 

M(t) m 



(12) 



(13) 



where Mi (m) is the mass of the i-th element released from a star 
with mass m. Because Ei(^ t) includes metals that already existed 
when the stars formed, these pre-existing metals must be subtracted 
in order to estimate the chemical yields, Pi(^ t). If the metallicity 
of stars with respect to the i-th element at their birth is , Ei can 
be divided into two terms, 



E^{^ t)=p,{^ t) + Z°E{^ t), 
where 

, dm 



P^{f^ t) 



y^{m)(|){m) 



(14) 



(15) 



M(t) 



and yi (m) indicates the mass of newly produced i-th element in a 
star with mass m. Portinari et al. (1998) show that the dependence 
of Mi and y; on initial stellar metallicity is small, so for simplic- 
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Figure 4. The cumulative SNe II rate. The lines types are the same as in 
Fig.H 



ity we neglect this dependence and use the values for solar initial 
metallicity throughout. 

In Figs.|4]and|5| we show Ru{^ t) and pi{^ t), respectively. 
In both cases, the solid and dashed lines denote the Kennicutt and 
top-heavy IMFs respectively. The initial time for which results are 
plotted corresponds to r (m„) and the final time to t{8Mq), where 
T{m) is the lifetime of a star of mass m. In Fig.|5| the lines indicate 
the yields of O, Si, Mg, Fe from top to bottom at the right hand 
side of the panel. In this paper, because instantaneous recycling is 
assumed for SNe II, we only use the values of chemical yields at 
t{8Mq). 

Note that we rescale the Mg yields from SNe II, pMg, 
upwards by a factor of 4, so that the Mg abundances pre- 
dicted by the model are more consistent with observed abun- 
dances in solar neighbourhood stars. This correction is stan- 
dard practice in chemical evolution modelling, and is jus- 
tified because the theoretical calculations of the yield of 
Mg are known to have a large (factor ^ 3) uncertainty 
jTinunesJ^oosk^^^^aveJ_^9^_jThom^^ 



[ 19981: IPortinari. Chiosi & Bressanlll998l: iLia. Portinari & Carrarc 
12002b . In reality, any correction to the yields yMg{m) is likely to 
depend on stellar mass, so the correction factor for the integrated 
yield pMg could be different for the two IMFs we use. In our model, 
the enrichment of galactic disks is dominated by stars formed with 
the Kennicutt IMF, while the a-element enrichment of the ICM 
is predicted to be dominated by stars formed with the top-heavy 
IMF, and in fact, comparison with solar neighbourhood abundances 
favours a lower correction factor (~2-3) for the Mg yield than do 
the ICM abundances. However, for simplicity we have assumed the 
same correction factor for both IMFs. 

The standard scenario (Whel an&Ib'en1ll973h for SNe la 
is that they occur in binary systems, in which a C-O white 
dwarf produced by evolution of the primary star accretes gas 
from the secondary star when the secondary overflows its Roche 
lobe at the end of its own evolution. This gas accretion drives 
the mass of the primary white dwarf over the Chandrasekhar 



o 
J 




Log[t/yrs] 



Figure 5. The cumulative chemical yields of SNe II, Pi(^ t). The solid 
and dashed lines denote the Kennicutt and top-heavy IMFs respectively. 
For each case, the lines indicate O, Si, Mg and Fe in sequence from top to 
bottom at t = t{8Mq) = 4 x 10'' yr. 



limit, causing it to explode. We model this using the scheme of 
Greggio & Renzini 1 1983), with parameters updated according to 
jPortinari. Chiosi & BressanI jl998l) . The progenitors of SNe la are 
assumed to be binary systems with initial masses in the range 
niB.iov < m.B < ms.up, where m_g = mi + m.2, and mi and 
m2 are the initial masses of the primary and secondary (m2 ^5 mi) 
respectively. We assume m_B,up = 2mi,max, where mi, max is 
the largest single-star mass for which the endpoint is a C-O white 
dwarf. The binary star systems which are the progenitors of SNe la 
are assumed to have an initial mass function A(j>{mB), where </>(m) 
is the same as for the single-star IMF. The distribution of mass frac- 
tions for the secondary, ^ — m2/m_g, is assumed to have the form 
(normalised over the range < pi < 1/2) 



-7K 



(0 < ^ < 1/2). 



(16) 



For a single generation of stars formed at t = 0, the cumulative 
number of SNe la explosions up to time t is then given by 



(jnB) 



1/2 



a(t) 



dmB 
rriB 



where the lower limit 

niB 



Mmin(i) = : 



ms - ms,up/2 



rriB 



(17) 



(18) 



is set by the conditions that the secondary has evolved off the main 
sequence and that mi = ttib ~ m2 ^ mi, max- Following Porti- 
nari et al. (1998), we adopt ms,iow = 3Mq, ms, up = 12Mq 
and 7 = 2. The factor A is usually chosen (for given tubaovj and 
7) to reproduce the observed ratio of rates of SNe la to SNe II in 
spiral galaxies. Following Portinari et al., we adopt A = 0.07. This 
observational normalisation is derived for an IMF similar to Ken- 
nicutt, but for simplicity we assume that the same normalisation 
applies to the top-heavy IMF also. 

In Fig.|6| we show the cumulative number of SNe la as a func- 
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Figure 6. The cumulative SNe la rate, i?ia- The hue types are the same 
as those used in Fig. |2] Note that the SNe la rate for the top-heavy IMF 
is smaller than that for the Kennicutt IMF, because the top-heavy DVIF has 
fewer intermediate mass stars. 

tion of time, Ria{^ t), for the Kennicutt IMF (solid line) and the 
top-heavy IMF (dashed line). The number of intermediate mass 
stars responsible for SNe la is larger in the case of the Kennicutt 
IMF than for the top-heavy IMF (as shown by Fig. |2j. This is in 
contrast to the situation for SNe II. The yield of metals from SNe la 
is computed as 

Pla,»(sS t) = M^i?Ia« t), (19) 

where M^'' is taken from the W7 model of iNomoto et alJ <1997h 
and takes the values ofO.143, 8.5x 10"^, 0.154 and 0.626 Mg for 
O, Mg, Si and Fe respectively. 

Note that if the IMF includes a fraction of brown dwarfs (with 
T > 1), the values in equations (TT}, GU, (BJ, {n} and (19) must 
be multiplied by 1/T, because T > 1 corresponds to having a 
fraction 1/T of visible stars with md < m < m„. 

3.2 The metal enrichment of cold and hot gas 

We treat the metal enrichment due to SNe II in the same way 
as CLBF, except that the instantaneously returned mass fraction 
and yields (denoted R and p in CLBF) are replaced here by 
r[8M0]) and Pi{^ t[8A'Iq]). The mass released from evolved 
stars after t{8Mq) and the metals produced by SNe la are com- 
puted at every timestep. The gas and metals recycled from disk 
stars are put into the cold gas in the galaxy, while those from bulge 
stars are put into the hot gas in the halo. This picture is consistent 
with many evolutionary models of disk galaxies, in which metals 
from SNe la are put into the gaseous disk; observations suggest 
that elliptical galaxies do not keep the metals released from SNe la, 
because the metal abundances in the interstellar medium of ellip- 
ticals are either smaller than or at best si milar to solar abundances 
jMatsushita. Ohashi & Makishimal2000h . 

The mass released from stars during the time interval ti <t < 
ti+i is computed from 



AM= [ [E{^ U+i - t') - U - t')]ip{t')dt' , (20) 
Jo 

where is the star formation rate at time t for the galaxy under 
consideration. In order to solve this equation, the star formation his- 
tories in both modes of star formation, quiescent and starburst, must 
be recorded separately for individual galaxies. Due to the memory 
overhead involved, the star formation histories are rebinned into 30 
linear steps in time. We have confirmed that our results are insensi- 
tive to the precise number of steps used. 

As well as the returned mass, the metals produced by evolving 
stars, including SNe la, are computed in a similar way. The mass 
of the j**^ element released from evolved stars and SNe la over the 
time interval [ti,ti+i} is 

AM,= Ml^lRi^s^ U+i - t') - Ri4^ U - t')]^{t')dt' 

Jo 

+ r'^'zit')[E{^ ti+1 ~ t') - U - t')]i];{t')dt' , (21) 
Jo 

where Z(t) is the metallicity of stars at their birth at time, t; this 
quantity is stored in the same fashion as the star formation histories. 
The metals are also assumed to be put into cold and hot phases 
along with the recycled gas. 



4 THE GALAXY LUMINOSITY FUNCTION 

In this paper we consider three principal models: the fiducial model 
of CLBF (which assumed — 0.02) and two new models which 
have fib — 0.04. The first of these new models is super-wind model, 
in which gas ejection from the halo is invoked to suppress the for- 
mation of bright galaxies. In the second new model, which we refer 
to as the conduction model, cooling in massive halos is suppressed 
by thermal conduction. In the CLBF model, all stars form with a 
Kennicutt IMF. In the siiperwind and conduction models, the Ken- 
nicutt IMF is adopted for quiescent star formation in disks and a 
flat IMF is used in starbursts. The superwind model is similar to 
the model used by Baugh et al. (2004b). We also consider a variant 
of the superwind model (which we denote superwind/Kennicutt) in 
which all star formation produces stars with a Kennicutt IMF. We 
list the parameters for the different models in TableQ 

An important constraint on the models is that they should be 
consistent with the observed present-day galaxy luminosity func- 
tion. In Fig. we compare the model predictions for the present- 
day A'-band luminosity function with observational estimates. The 
predictions of the superwind model with Kccap = 600 km s~^ 
are shown by the solid line, those of the conduction model with 
Kond = 100 km s^^ by the short dashed line, and the luminosity 
function of the CLBF model is plotted using a long dashed line. 
The symbols show observational data. Although it is challenging 
for models with the currently favoured value of fib to match the ob- 
served luminosity function as well as in the CLBF model, as shown 
by Benson et al. ( 2003a), the superwind and conduction models are 
in broad agreement with the observations. 

In the conduction model, the value of the parameter V^^-^^ 
which is required to produce a good match with the observed lu- 
minosity function seems unphysically small. The value of l^ond = 
lOOkms^^ requires a thermal conductivity that is many times 
larger than the Spitzer value. As explained in Section 2, processes 
other than thermal conduction could also contribute to the suppres- 
sion of cooling, so we retain the freedom to adopt such a low value 
for this parameter. 
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Table 1. Model parameters. Column 1 gives the name of the model, column 2 gives the baryon density parameter, column 3 gives F, the CDM power spectrum 
shape parameter Column 4 gives the dark matter halo mass function used; J and PS denote mass functions given by Jenkins et al. 1 2001 ) and Press & SchechteJ 
1 19741) respectively. Columns 5 and 6 give the parameters that specify the quiescent star formation timescale, as defined by eqn. Q. Note that in the CLBF 
model, the star formation timescale scales with the dynamical time of the disk, Tjigij. Column 7 gives the value of the V^ot parameter that determines the 
strength of the standai'd SNe feedback (eqn.|3). Columns 8 and 9 describe the superwind feedback; column 8 gives the parameter that sets the strength of the 
superwind (eqn.|8) and column 9 gives the threshold circular velocity of halos above which gas ejected by a superwind can be recaptured. Column 10 gives 
the circular velocity of halos above which cooling is suppressed in the conduction model. Column 1 1 gives the IMF assumed in bursts; note that all models 
assume a Kennicutt IMF for quiescent star formation. Finally, column 12 gives the ratio of total stellar mass to the mass produced in luminous stars. 
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Figure 7. The K-band galaxy luminosity function at 2: = 0. The model 
predictions are shown by the lines as indi cated by the key. Th e symbols with 
error bars show ob servational data_from lColeet alj l200ll filled circles), 
lHuangetal]<20O3l stars). and iBell et alj<2003l open circles). 
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Figure 8. The K-band luminosity function predicted in the superwind 
model for different values of Viocap, the threshold halo circular velocity 
for the recapture of gas previously ejected by a superwind (as indicated by 
the key). The fiducial superwind model has Vrecap = 600kms~^ and is 
shown by the solid line. The observational data are plotted using the same 
symbols as in Fig. 171 



In Fig. |8|we show how the luminosity function in the super- 
wind model depends on the parameter Kccap, the threshold circu- 
lar velocity above which halos recapture gas previously ejected in 
a superwind. Clearly if Vrocap is too low, then too many very bright 
galaxies are produced. On the other hand, the galaxy luminosity 
function becomes insensitive to increases of this parameter above 
Kccap ~ 600 km s~^. However, as we will show later, the value 
of Kccap has a strong impact on metal abundances in the ICM. 
The predicted ICM metal abundances are too low unless the metals 
ejected in the superwind can be recaptured by cluster-sized halos. 
Thus, we require a value Kccap — 600 km s~^ in order to repro- 
duce both the observed field galaxy luminosity function and the 
metal abundance of the ICM. 



5 PRESENT-DAY INTRACLUSTER MEDIUM 

5.1 Gas mass-to-stellar luminosity ratio and hot gas fraction 

A basic test for any model of the chemical enrichment of the ICM 
is that it should reproduce the observed ratios of the ICM mass to 
total stellar luminosity and to total mass in clusters. This is because 
the metals are produced in stars, so a model which produces the 
observed ICM metallicity is invalid unless it predicts the correct 
number of stars at the same time. 

In Fig.|9| we plot the ratio A/hot /^bj of the mass of the ICM, 
Mhot, to the total fej-band stellar luminosity of the galaxies, L^j, 
as a function of the ICM temperature T. We show predictions for 
the superwind, conduction and CLBF models. In the models, the 
ICM temperature is assumed to equal the halo virial temperature, 
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J-vir = -TT— Vr. ~ O.lkeV 



2 fef 



V lOOOkms" 



(22) 



\ (a) superwind 



while for the observations we use the measured X-ray tem- 
perature. The relationship between cluster mass and virial 
temperature in the ACDM model at 2: = is then 
T^ir = 6.4keV(M/10"fe~^Mffi )^/^ We compa r e the mod- 
els with observational dat a from ISanderson et alj i2003h and 
ISanderson & PonmanI i2003h . We have combined the data from in- 
dividual clusters in bins in X-ray temperature, and computed the 
mean and dispersion in each bin, weighted according to the errors 
on the individual datapoints. The error bars on the plotted data- 
points are calculated from the dispersion in each bin, and do not 
include any other measurement errors. 

The conduction and CLBF models predict a ratio A/hot /it, 
which increases only slightly with cluster temperature. The super- 
wind model predicts a somewhat larger increase, because of the 
expulsion of gas from lower-mass halos, but again is almost flat 
for halos with Vc > Kecap, corresponding to T ~ l.lkeV which 
is indicated by the vertical dashed line in the top panel. The mod- 
els are in approximate agreement with the observational data for 
T > 2keV, but for lower temperatures, the observations suggest 
ICM fractions significantly smaller than what the models predict. 
However, there are still important uncertainties in the observational 
data (not included in the error bars plotted) resulting from the need 
to extrapolate X-ray measurements out to the cluster virial radius, 
so the ICM fractions in cooler clusters are not yet securely deter- 
mined as shown below. 

We note that the CLBF model predicts very similar A/hot /^tj 
values to the superwind and conduction models, even though it has 
a global baryon fraction only half that of the other models. This 
agreement comes about because the CLBF model predicts a lumi- 
nosity density which is a factor ~ 2 lower in all environments com- 
pared to the other models. 

In Fig. llOl we show the ratio of the ICM mass to the total clus- 
ter mass including dark matter as a function of cluster temperature. 
The dots and crosses indicate the hot gas fractions for the super- 
wind and the conduction models, respectively. The dashed hori- 
zontal line corresponds to the cosmic baryon-to-total mass ratio; 
the model points would all lie on this line if all baryons were re- 
tained in the cluster in the hot gas component. We also show obser- 
vational data from Sanderson et al. (2003), binned in temperature 
as in Fig.|9| The superwind model reproduces the qualitative trend 
seen in the observations that the gas fraction drops in low tempera- 
ture clusters. However, we caution that the observational estimates 
rely on extrapolating both the ICM and dynamical masses from the 
outermost measured X-ray datapoint out to the halo virial radius. 
Future more accurate X-ray observations will allow more precise 
tests of models for gas ejection from halos. 



5.2 Metal abundances 

We now compare the ICM metallicities predicted by our mod- 
els with observational data. Our model predicts global metallic- 
ities, averaged over the entire ICM within the cluster virial ra- 
dius. In contrast, the ICM metallicities measured from X-ray data 
are often only representative of the central regions of the clus- 
ter, because these dominate the X-ray emission. Some clusters, 
known as cooling flow clusters, show metallicity gradients for 
some elements, in the sense that metallicities increase towards 
the cluster centre. It is therefore essential to correct the obser- 
vational data to estimate mean metallicities averaged over the 
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Figure 9. The ratio of ICM mass to bj-band luminosity of stars contained 
in a cluster, Mjiot/^bj, plotted as function of the hot gas temperature. 
The dots show the predictions of: (a) the superwind model, (b) the con- 
duction model, and (c) the CLBF model. The crosses in all panels denote 
the binned averages of M\, r,t/Lh-, estimated using observ ational data from 
ISanderson et alj <2003ft andjsanderson & Ponman ('2003). The dashed ver- 
tical line in panel (a) shows the temperature corresponding to a halo circular 
velocity V^ocap = 600 km s~^, above which ejected gas is recaptured. 



entire cluster, before comparing with the models. We compare 
our model predic tions with the followin g observational data on 
ICM abundances : De Grandi et al. (2004, BeppoSAX) for Fe [see 
also Irwin & Bregman t2001) and De Grandi & Molendi l2Q02) 1; 
IPeterson et alJ <2q03l X M M-Newton) for O Fe, Mg and Si [see 
also iTamura et al.' ^004*)! : 'Fukazawa et al J j 19981 ASCA) for Fe 
and Si; and Baumgartner et al. ( 2003, ASCA) for Si. We explain in 
the Appendix why we have omitted some observational data from 
our compari son. We use the results of the systematic study of iron 
gradients bv lDe Grandi et aO i2004b to convert the measured cen- 
tral metallicities of iron and silicon to global average values, as- 
suming that silicon has the same gradient as iron, as is suggested 
by observational data. We present the details of how we make these 
corrections in the Appendix. The size of these estimated correc- 
tions is up to a factor of two. Oxygen and magn esium do not show 
metallicity gradients framura etal 200l"2004f), so we assume that 
the global abundances of these elements are equal to the measured 
central values. We note that previous studies (e.g. de Lucia et al. 
2004) have not included any correction for metallicity gradients 
when comparing models with observed ICM abundances. 



We first show how the predicted metal abundances in the ICM 
depend on the IMF adopted for starbursts. Fig. II II shows the O, 
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Figure 11. The metal abundances of the ICM. Each panel shows a different element: (a) [O/H] (b) [Fe/H] (c) [Mg/H] (d) [Si/H]. Abundances are normalized 
to solar values, and shown on a logarithmic scale. The predictions of the superwind model are shown by dots in the standard case where starbursts have a 
top-heavy IMF, and by crosses for the variant model in which all star formation takes place with a Kennicutt IMF. The open triangles in each panel denote 
observational data from XMM-Newton i Peterson et al. 2003), and the open squares data from ASCA ( Fukazawa et al. 1998); in panel (b), the open circles show 
data from BeppoSAX ( De Grandi et al. 2004), while in panel (d), the open circles show A5CA data from Baumgartner et al. 1 2003). All of the datapoints are 
measurements for individual clusters, apart from Bau mgartner et aU 1.20031) . which are average values in bins of temperature. The observational data have been 
corrected for abundance gradients as described in the text. 



Fe, Mg and Si abundances (relative to solar values) as functions 
of cluster temperature for the superwind model. The dots show 
the predicted values for the standard superwind model with a top- 
heavy IMF (x = 0) in bursts and a Kennicutt IMF in disks, and the 
crosses show a variant superwind model in which disks and bursts 
both have a Kennicutt IMF. For clusters with T > IkeV, the abun- 
dances of the a-elements O, Mg and Si are 3-4 times higher in the 
case of the top-heavy IMF in bursts, while the abundance of Fe is 
only about 60% higher. This difference reflects that fact that the a- 
elements are produced mostly by SNe II, and the number of SNe II 
is larger for a top-heavy IMF, while Fe is produced by both SNe la 
and SNe II, and the number of SNe la is actually reduced for the 
top-heavy IMF. In fact, the models predict that most Fe is produced 
by SNe la for a Kennicutt IMF, but by SNe II for the top-heavy 
IMF 

When in Fig.fTTI we examine the dependence of ICM abun- 
dances on cluster temperature (or mass) in the superwind model, we 
see that the abundances of the different elements all show a break 
at T ~ 1.1 keV. This temperature corresponds exactly to the halo 
circular velocity Viccap = 600 km s~^ above which gas and met- 
als ejected by superwinds are recaptured. Halos with Vc > 'Kccap 
effectively trap all of the metals ever produced by the galaxies they 



contain. For clusters hotter than the break, the abundances are pre- 
dicted to be almost constant. The jump in abundances across the 
break is about a factor 2 for O, Mg and Si, but is much smaller 
for Fe, especially if a Kennicutt IMF is used in bursts. This again 
reflects the fractions of these different elements produced in either 
SNe II or SNe la. Since SNe II also drive the superwinds, while 
SNe la are typically delayed until after the strong superwind phase 
has subsided, superwinds are much more effective in expelling 
from halos the heavy elements produced in SNe II than those pro- 
duced in SNe la. If Kccap is varied, then the metal abundances in 
clusters with Vc > Vrocap are insensitive to the value of Kocap 
for Kccap 400 km s^^, but increase slightly if Kccap ^ 400 
km s^^, because more stars are produced in that case (c.f. Fig|8j. 

Fig. ll ll also shows the observed ICM abundances of O, Mg, 
Si and Fe, corrected to global average values as we have described. 
We see that the superwind model with a top-heavy IMF in bursts 
agrees well with the observed abundances of O, Mg and Si, while 
the superwind model with a Kennicutt IMF predicts abundances of 
these elements which are too low by factors 2-3. For Fe, the model 
with the top-heavy IMF again agrees better with the observations, 
although the abundance is still slightly low. We note that the cur- 
rent abundance measurements for O, Mg and Si are all for clusters 
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Figure 10. The ratio of hot gas mass to total mass, including dark matter, 
as a function of the hot gas temperature. The dots and crosses show the 
predictions of the superwind model and the conduction model, respectively. 
The point s with error bars show the binned averages of the observational 
data from ISanderson et alj l2003l) (we use their estimates of gas fraction 
wihthin the virial radius). The horizontal dashed line indicates the universal 
ratio of baryon to total mass density estimated from Qo and Sl^ . The vertical 
dashed line shows the temperature corresponding to a halo circular velocity 
V, 

recaptured. 



recap = 600 km s ^ , above which ejected gas in the superwind model is 



hotter than the break predicted in the models. However, the value of 
Kecap cannot be significantly increased above our standard value of 
600 km s~^ without moving the abundance break into the tempera- 
ture range where the observations indicate a plateau in abundances, 
thus conflicting with the data. The measured Fe abundances, which 
extend down to somewhat lower cluster temperatures, do in fact 
suggest a break close to the temperature predicted by the super- 
wind model, but more measurements for cool (T < IkeV) clusters 
are needed to confirm this. It is clear that more abundance mea- 
surements in cool clusters, especially of O, Mg and Si, would put 
stringent constraints on our model of ejection of gas and metals 
from halos by superwinds. The models also predict a very small 
scatter in global ICM abundances at a given cluster temperature, 
especially above the break. The observed abundances show a larger 
scatter, but most of this results from measurement errors. 

In Fig. 1 121 we show the ICM metal abundances predicted by 
our other two models: the conduction model with a top-heavy IMF 
in bursts, and the CLBF model with a Kennicutt IMF for all modes 
of star formation. The models are compared with the same observa- 
tional data as in Fig. ll II For both of these models, halos retain all 
of their metals, and consequently the predicted ICM abundances 
are almost constant with cluster temperature. The CLBF model 
predicts abundances too low compared to observed values, by fac- 
tors 2-3 for O, Mg and Si, and by around 40% for Fe, and is thus 
ruled out. The conduction model predicts abundances in much bet- 
ter agreement with observations, although slightly too high for O 
and Si. The abundances in the conduction model are slightly higher 
overall than those in the superwind model. 



It should be remembered that we have applied a correction to 
the yield of Mg from SNe II predicted by stellar models, in order 
to better match observational constraints. The model yields have 
been multiplied by a factor of four throughout, a value which is a 
compromise between that required to reproduce Mg abundances in 
galactic disks (specifically, in the solar neighbourhood) and in the 
ICM; a separate correction factor could be applied to each IMF be- 
cause of the mass-dependent uncertainties in estimating theoretical 
yields. No correction factors have been applied to the yields of the 
other elements: O, Fe and Si. 



5.3 Abundance ratios 

Examining the abundance ratios between different elements in the 
ICM provides an alternative method to constrain the form of the 
IMF, in which the effects of uncertainties in the ICM to stellar 
mass ratio cancel out. In Fig. II3I we show the abundance ratios 
of different q elements to iron, [O/Fe], [Mg/Fe] and [Si/Fe], for 
the standard superwind model, with a top-heavy IMF in bursts, and 
the variant superwind model, wit h a Kennicutt IMF in bursts. We 
also show observational data from lPeterson et all i2003h . which in 
the case of Si and Fe have been corrected for abundance gradients 
in the way already described. We see that the predicted O/Fe and 
Mg/Fe ratios agree better with observations for the model with the 
top-heavy IMF, while on the other hand the Si/Fe ratio agrees better 
with the model having a Kennicutt IMF in bursts. The overall level 
of agreement of the abundance ratios with observational data in the 
top-heavy IMF model could be improved if the Fe abundances in 
the model were slightly larger (this can also be seen in Fig. ll 1> . or 
if the observed Si abundances were slightly larger. The latter would 
occur if the abundance gradients in Si were weaker than we have 
assumed, since we would then infer larger global Si abundances 
from the same meas ured central values. However, recent data from 
Tam ura et alj i2004) . in which metal abundances are measured out 
to r < 500/i~^ kpc from cluster centres, seem to support our orig- 
inal assumption that the radial gradients in Si and Fe are the same. 

Fig. 1141 is the equivalent of Fig 1131 for the conduction and 
CLBF models. Since the conduction model predicts slightly higher 
Fe abundances than the superwind model, the agreement with the 
observed abundance ratios improves. The CLBF model fails to re- 
produce the observed abundance ratios. The models predict very lit- 
tle variation in the abundance ratios with cluster temperature, apart 
from the break in the superwind model where recapture sets in. 
They also predict very small scatter for the abundance ratios at a 
given cluster temperature. 



5.4 Metallicities of stars and gas in clusters 

In Fig. [B] we compare the predicted mass-weighted mean iron 
metallicities for the different baryonic components in clusters: stars 
in galaxies, cold gas in galaxies, and hot intracluster gas. Results 
are shown for the standard superwind and conduction models as 
functions of cluster temperature. In both models, the mean stellar 
and cold gas metallicities are close to solar, while the hot gas metal- 
licity is 2-3 times lower. The mean total metallicity is only slightly 
larger than the hot gas metallicity, reflecting the fact that most of 
the baryons are in the hot gas component. For the superwind model, 
the onset of recapture of gas by halos produces the expected jump 
in the hot gas and total metallicities, but does not produce any cor- 
responding feature in the metallicities of the stars or cold gas. This 
is because the enrichment of these components has almost finished 
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Figure 12. The metal abundances in the ICM for the conduction model (dots) and the CLBF model (crosses). The observational data are plotted using the 
same symbols as in Fig. II II 



by the formation epoch of the cluster and the epoch of recapture of 
expelled metals (see Section 7). 

Fig. 116! shows the fractions of the total mass of iron which 
are in stars, cold gas, hot gas and 'lost' gas which is expelled from 
halos by superwinds, for the same models as in Fig.[T5l In all cases, 
more of the iron is in the hot gas than in galaxies. In the superwind 
model, superwinds put about 40% of the iron outside halos for low 
mass clusters below the threshold for recapture. Thus, according to 
this model, a significant mass of metals should be observed in the 
surroundings of low temperature clusters. 

Thus, in summary, from the ICM results at 2: = 0, we con- 
clude that a top-heavy IMF with a: ~ in starbursts is required 
for the model predictions to match the observed ICM metallic- 
ities, as long as we adopt the Kennicutt or a similar IMF for 
disk star formation. The use of a Kennicutt or a similar IMF for 
disk star formation is supported by many observations and mod- 
els of^ol^jieighbourhood stars, i ncluding a semi-analytic model 
by iNaeashima & Okamotd i2004) in which a Salpeter-like IMF 
provides good agreement with the observed metallicities of solar- 
neighbourhood stars. We have also confirmed that while using 
the Arimoto-Yoshii IMF ( Arimoto & Yoshii 1987) with x — 1 
for both star formation modes substantially increases metallicities 
compared with the Kennicutt IMF, a combination of the Kennicutt 
IMF for quiescent star formation and the Arimoto-Yoshii IMF for 
starbursts does not reach the observed metallicities. 



6 ICM ABUNDANCES IN HIGH-REDSHIFT CLUSTERS 
6.1 Evolution of [Fe/H] 

With the advent of the new generation of X-ray satellites, measure- 
ments of ICM metallicities are now becoming possible in high- 
redshift clusters, thus allowing direct tests of ICM evolution. In 
Fig. ll7l we show what our models predict for the evolution of the 
iron abundance in the ICM with redshift in the range < z < 1.3. 
At each redshift, we calculate the mean Fe abundance for clus- 
ters in two different ranges of temperature, 3 < T < 5keV and 
5 < T < 8.5keV. We show results for the standard conduction 
and superwind models. In the sam e figure we also show observa- 
tional data from lTozzi et al .'l2Q03). based on measurements of the 
iron abundances in 18 distant clusters with 0.3 z < 1.3 using 
Chandra and XMM-Newton. We have binned this data in redshift 
and in the same ranges of temperature as used for the models. Note 
that we have no information about metallicity gradients for these 
distant clusters, so we have not applied any corrections to the mea- 
sured abundances to get global values. We also show observational 
data for z = Q clusters binned in the same ranges of temperature. 
The z = data are corrected for metallicity gradients as previously 
described. 

Our models predict a weak decline in the Fe abundance with 
increasing redshift, by about 30% over the range < 2; < 1.3, and 
furthermore this decline is pre dicted to be the sa me for both low 
and high temperature clusters. ItozzI et all i2003h concluded from 
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Figure 13. The abundance ratios of a-elements to iron. Each panel corre- 
sponds to a different ratio: (a) [O/Fe]. (b) [Mg/Fe]. (c) [Si/Fe]. All ratios are 
relative to solar values. The predictions of the superwind model are plotted 
using dots for the standard case in which starbursts have a top-heavy IMF, 
and crosses when all stai' formation is with the Kennicutt I MF. The open 
triangles show observational data from lPeterson et al] l2003l) . corrected for 
radial gradients. 



their own data that the evolution in Fe abundance with redshift is 
negligible, but we see that the data for the hotter (5 < T < 8.5keV) 
clusters are also quite consistent with the weak decline predicted 
by our models. The data for the cooler (3 < T < 5keV) clusters 
suggest higher Fe abundances and a weak increase with redshift 
relative to the hotter clusters, both of which would conflict with 
our model, but the uncertainties are very large. 



6.2 Evolution of [O/Fe] 

The abundances of elements other than Fe have not yet been mea- 
sured in high-redshift clusters, so we simply present predictions 
for the evolution of O abundances from our model. For clusters se- 
lected in the same temperature range at different redshifts, the mean 
O abundance is predicted to decline with redshift even more weakly 
than for Fe, by only ^ 10% in over the range < z < 1.3. The 
difference in behaviour is caused by the difference in the timescales 
for enrichment by SNe la and SNe II. The O/Fe ratio in the ICM 
is therefore predicted to increase with redshift. This is shown in 
Fie. 1181 For the superwind model, there is a significant difference 
in O/Fe between low (0.7-1 keV) and high (3-8.5 keV) temperature 
clusters, caused by recapture of superwind ejecta in clusters with 
T > l.lkeV. 
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Figure 14. The abundance ratios in the conduction model (dots) and the 
CLBF model (crosses). The data are the same as plotted in Fig. 1131 
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Figure 15. The iron metallicities, [Fe/H], of stars, cold gas and hot gas, 
shown by the solid, dotted and short-dashed lines, respectively. The mass- 
weighted, mean total metallicities are shown by the long-dashed line. Panel 
(a) shows results for the superwind model and (b) for the conduction model. 
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Figure 19. The evolution of the main progenitors of present-day clusters, as defined in the text. The left, middle and right panels in each row show the evolution 
of clusters with low (0.7-1 keV), intermediate (1.5-4 keV) and high (4-8 keV) present-day temperatures respectively. Top row: The evolution of halo mass of 
the main progenitor, relative to the present-day cluster mass. Second and third rows: The evolution of iron abundance [Fe/H] in the ICM of the main progenitor, 
in the superwind and conduction models respectively. Fourth row: The evolution of the oxygen abundance [O/H] in the ICM of the main progenitor, in the 
superwind model. In all panels, the line style changes from dotted to solid once the formation redshift of the cluster is passed. 
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Figure 17. The redshift evolution of [Fe/H] in the ICM. The lines show 
the mean abundances predicted by the models for two different temperature 
ranges, 3-5keV (dotted) and 5-8.5keV (solid), for both the superwind and 
conduction models (lower and upper pairs of lines, as labelled). The eiTor 
bars on the hnes show the Icr scatter. The symbols show the observal data, 
with means and Icr scatter comp uted in each bin. P entagons indicate the 
data for high redshift clusters from lTozzi et alJ i20O3h . circ les, triangles and 
squai'es are mean metallicities of clusters at z = given b y lPe Grandi et alj 
( 2004), Peterson et al. ( 2003) and Fukaz awa et alJHQQSh . respectively. The 
symbols for z = are given small offsets for clarity. The observational data 
are binned into the same temperature ranges as the models. 
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Figure 18. The redshift evolution of [O/Fe] in the ICM. The lines show 
the predictions for the conduction and superwind models respectively. In 
the case of the superwind model, the sample is divided into two tempera- 
ture ranges: high (3-8.5 keV; solid line) and low (0.7-1 keV; dotted line) 
temperatures. 



7 METAL ENRICHMENT HISTORIES 

In this section, we investigate what our model predicts for the 
enrichment histories of individual clusters. For this purpose, we 
divide clusters into three populations based on their present-day 
ICM temperatures: low (0.7-1 keV), intermediate (1.5-4 keV) and 
high (4-8 keV), corresponding to cluster masses 3.6 — 6.2 x 10^'^, 
1.1 - 5.0 X 10" and 0.5 - 1.4 x lO^^/i~^Af0 respectively. These 
ranges are chosen to illustrate different behaviours regarding recap- 
ture of gas by halos in the superwind model. For the low temper- 
ature clusters, a significant quantity of metals and gas have been 
expelled from the progenitors of these halos without recapture. The 
intermediate temperature clusters experience recapture of gas and 
metals at low redshifts, so we can see how the recapture process 
affects the evolution of the ICM metallicity. The high temperature 
clusters have already recaptured gas and metals at high redshifts. 

In Fig.[l9l we show the evolution of the halo mass and ICM 
metallicity for the main progenitors of present-day clusters, for a 
randomly chosen set of clusters in each range of present-day tem- 
perature. We define the main progenitor of a present-day cluster at 
any earlier epoch as follows (Lacev & Cole 1993): starting from 
the present-day halo, we follow the halo merger history back in 
time, and at each merger event, we follow the branch in the merger 
tree corresponding to the more massive progenitor involved in that 
merger event. We define the formation redshift 2/ of a cluster as the 
earliest redshift at which the main progenitor has at least half of the 
present cluster halo mass liLacev&ColeHl993l) . For z ^ z/, the 
main progenitor is also the most massive progenitor, but at earlier 
times this may not be the case. 

The top row of panels in Fig. 1191 shows the evolution of the 
halo mass of the main progenitor relative to the present-day halo 
mass, illustrating the spread in cluster assembly histories. We plot 
each cluster trajectory in Fig.^|as a solid line for z < Zf and as 
a dotted line for z > z/. In the ACDM model, clusters typically 
form around z ~ 0.5, but with a wide dispersion, and with more 
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massive clusters typically forming more recently. The medians and 
the 10% and 90% percentiles of the distribution of z/ values are 
(0.64,0.38,0.95) for low temperature, (0.57,0.27,0.79) for medium 
temperature, and (0.46,0.21,0.75) for high temperature clusters re- 
spectively. 

The lower panels of Fig. 1191 show the evolution of the ICM 
abundances in the main progenitor of each cluster, for the same set 
of cluster halo assembly histories as in the top row. The second and 
third rows show the Fe abundance in the superwind and conduction 
models, while the fourth row shows the O abundance in the super- 
wind model only. The Fe abundances typically evolve only weakly 
with redshift over the range < z < 1.5, which is similar to the 
behaviour seen in Fig. 1171 The dispersion in the final metallicity 
and in the metallicity evolution is larger in the lower temperature 
clusters. For the intermediate temperature clusters, the dispersion 
in the metallicity evolution is larger in the superwind model than in 
the conduction model, due to the ejection and then recapture of gas 
in the former case. These recapture events can be seen as abrupt 
upward jumps in metallicity in panel (e). The evolution of the O 
abundance in the superwind model shows a somewhat larger dis- 
persion than the evolution in Fe abundance. The jumps in O metal- 
licity due to recapture of superwind material are larger than for Fe 
(compare panels (k) and (e)), because the superwinds are preferen- 
tially enriched in products of SNe II rather than SNe la. Some of 
the evolution tracks for the O abundance show a decline in mean 
ICM abundance with increasing time. This behaviour can be un- 
derstood as follows: the dispersion in O abundance between halos 
of the same mass increases with decreasing halo mass, due to vari- 
ations in the number of ~ L, galaxies within a halo resulting from 
the stochastic nature of the halo assembly and galaxy formation 
processes. As the halo mass builds up, the metallicity tends to con- 
verge towards the average value for clusters with that present-day 
temperature: progenitor halos with above-average ICM metallicity 
tend to merge with other halos having lower metallicity, and the 
resulting mixing of the hot gas components produces the declining 
tracks seen for instance in panel (j). 

In Fig.|20l we show the build-up of the total mass of metals 
found in the present-day ICM for the same set of clusters as in 
Fig. II9I We calculate these evolutionary tracks by summing over 
the ICM metal masses of all of the progenitors of a present-day 
cluster. We include in this sum metals which have been ejected by 
superwinds, if these metals have been recaptured by the present 
day. These tracks thus allow us to determine the redshifts at which 
metals in the ICM of the present-day cluster were actually produced 
(in the sense of being ejected from stars in supernova explosions). 
The top two rows of panels show the evolution of the iron mass 
in the superwind and conduction models, for different ranges of 
present-day cluster temperature. The bottom row shows the evo- 
lution of the oxygen mass in the superwind model. Most of the 
metals are typically produced well before the clusters are assem- 
bled, with the oxygen being produced even earlier than the iron. 
For example, in the superwind model, half of the O is produced by 
redshifts zq = (1-8, 2.1, 2.1) for the low, medium and high tem- 
perature clusters respectively, compared to zpe ~ (1-0, 1.3, 1.3) 
for producing half of the Fe. The corresponding redshifts in the 
conduction model are not much different: zo ~ (1.7, 1.8, 1.9) and 
ZFe = (1.1, 1.2, 1.2). Thus, the ICM enrichment occured mostly 
in many smaller progenitors of a present-day cluster. We also see 
that metal production typically occured slightly earlier in high mass 
than low mass clusters, even though high mass clusters were assem- 
bled (in the sense of reaching half of their present mass) slightly 
later. 



We can compare our results for the redshift of metal produc- 
tion in clusters with the results from t he semi-analytical model of 
iDe Lucia. Kauffmann & White 1 12004). De Lucia et al. use the in- 
stantaneous recycling approximation throughout, so their results 
can be considered to apply to oxygen, but are not valid for iron, 
for which the delayed production by SNe la is important. Further- 
more, they only compute the ICM evolution for the single exam- 
ple of a lO^^/i~^M0 cluster. For three different feedback/ejection 
schemes, they calculate that 60-80%, 35-60%, and 20-45% of the 
metals are incorporated into the ICM at redshifts larger than 1 , 2, 
and 3, respectively. Note that in their definition, metals are "incor- 
porated" into the ICM either when they are ejected from galaxies 
directly into the ICM, or when metals ejected from the cluster are 
recaptured. Comparing with the predictions for high temperature 
clusters in our own models, we find that the fraction of oxygen in- 
corporated into the ICM at redshifts larger than 1, 2 and 3 is 47%, 
23%, and 10% respectively for the superwind model, and 79%, 
51%, and 24% for the conduction model. Thus in our superwind 
model, O is incorporated into the ICM later than in the models of 
De Lucia et al., while in our conduction model, the incorporation 
happens at similar redshifts. (Note that the "metal incorporation" 
redshifts given here are different from the "metal production" red- 
shifts derived from Fig. |20| since the latter include ejected met- 
als among the metals produced by a given redshift, provided these 
ejected metals have been recaptured by the present day.) Since we 
assume the same ACDM cosmology as De Lucia et al., the differ- 
ences must arise from using somewhat different assumptions about 
gas cooling, star formation, feedback and gas ejection. 

Mergers of halos are expected to lead to extensive mixing of 
the ICM, which will tend to erase any metallicity gradients. The 
prediction of the models that typically most of the iron in the ICM 
is produced well before half of the cluster mass has been assem- 
bled then seems to make it harder to understand the iron abundance 
gradients seen in "cooling flow" clusters. However, since the pro- 
duction of iron by SNe la is long-delayed, it may still be possible to 
build up iron abundance gradients in the central regions of clusters 
by ejection of iron from the central galaxy after most of the cluster 
mass has been assembled. This problem needs to be investigated 
using gas-dynamical simulations of cluster formation. 

We have also calculated the build-up of metals within cluster 
galaxies, counting metals in both stars and in the cold gas. When 
we compute trajectories for the evolution of the total mass in met- 
als in galaxies summed over all cluster progenitors, relative to the 
present metal mass, we find that there is very little dispersion be- 
tween different clusters of the same mass (much less than for the 
ICM metal mass shown in Fig. l20> . The trajectories are also very 
similar for different ranges of cluster temperature, and for the su- 
perwind and conduction models. For clusters in the temperature 
range 0.7 < T < 8keV, half of the O mass in cluster galaxies is 
produced by redshifts zo = 2.6 — 2.8 for the superwind model, 
and by zo = 2.7 — 2.9 in the conduction model. For the Fe mass in 
cluster galaxies, the corresponding redshifts are zpe = 1.7— 1.8 in 
the superwind model, and zpe ~ 1.9—2.0 in the conduction model. 
The metals in cluster galaxies are thus produced significantly ear- 
lier than the metals in the ICM. 

Finally, we examine which galaxies contribute the most to 
cluster metal enrichment. Fig. 1211 shows the total amount of iron 
which galaxies have produced (either retained within the galaxy, 
injected into the ICM, or ejected from the cluster) as a func- 
tion of their A'-band luminosities at 2: — 0. The contribution of 
each galaxy is calculated by summing the iron production over 
its formation history. The distributions over luminosity dF{^ 
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Figure 20. (a)-(f): The evolution of the total mass of metals in the ICM, summed over all progenitors, relative to the present-day metal mass. Metals ejected 
by superwinds ai'e included if these are recaptured by the present-day. The top two rows show the evolution of the iron mass in the superwind and conduction 
models respectively, while the bottom row shows the evolution of the oxygen mass in the superwind model only. The different columns show the results for 
different ranges of present-day cluster temperature in units of keV as indicated by the labels. 



Lif ) /din Lif are then normalized to unit area. Results are shown 
for both the superwind and conduction models, and also split ac- 
cording to present-day cluster temperature. The solid line in each 
panel shows the iron production averaged over halos of all masses 
(including halos much smaller than clusters). These distributions 
are quite similar to the distributions of total stellar mass, which is 
expected since the mass of iron produced is closely related to the 
mass of stars formed. We also find that the normalized distribu- 
tions for the mass of iron injected into the ICM are very similar to 
the distributions in Fig.|^for the total mass of iron produced. 

The distributions of iron production in Fig. l21l for the clusters 
show pronounced peaks at the high luminosity end. As the cluster 



mass or temperature increases, the position of the peak moves up 
in luminosity, but its amplitude decreases. This peak is produced 
by the central galaxies in the cluster halos, and reflects the distri- 
bution of total stellar mass (Benson et al. 2003b). The peak is most 
prominent for the lowest temperature clusters, because mergers are 
most effective at concentrating the stellar mass into a single cen- 
tral galaxy in such systems. Since the position of the peak due to 
central galaxies depends on halo mass, it is smeared out when we 
average over all halos, producing a smoother distribution (see the 
solid lines in Fig. l2U . 

We find that in the superwind model, half of the iron is pro- 
duced by galaxies with luminosities Lk brighter than 1.1 x 10^^, 
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Figure 21. The fractional contribution of galaxies to iron enrichment as a 
function of /C-band luminosity, Li(. Iron in all baryonic components (stars, 
cold gas, hot gas or ejected gas) is included, and the distributions are nor- 
malized to unit area. Panel (a) shows results for the superwind model and (b) 
for the conduction model. The long-dashed, dotted, and short-dashed lines 
indicate high, intermediate, and low temperature clusters, respectively. The 
solid lines show the distributions averaged over halos of all masses. 



4.2 X 10^" and 3.6 x 10^°h-^LQ for low, medium and high 
temperature clusters respectively, corresponding to galaxy stellar 
masses 5.4 x lO", 2.2 x 10^" and 1.9 x lO^'/i^^M© respectively. 
The median galaxy masses and luminosities for oxygen produc- 
tion are almost identical to those for iron. In the conduction model, 
the median galaxy masses and luminosities for metal production 
are somewhat higher: for iron production, the median K-band lu- 
minosities are 1.2 x lO", 6.2 x lO"' and 5.4 x 10"/i"^Lq, 
and the median stellar masses are 7.3 X 10^°, 3.5 X 10^° and 
3.1 X lO^°/i~'^M0 respectively, with very similar values for oxy- 
gen production. When we look at metal production averaged over 
halos of all masses, we find median luminosities Lk for produc- 
ing Fe or O of 4.2 x 10^° and 5.5 x IO^/i^^Lq respectively for 
the superwind and conduction models, and median stellar masses 
of 1.8 X 10^° and 2.5 x 10^°/i"^Mq respectively. We thus see 
that the median K-band luminosity (or stellar mass) for metal pro- 
duction is similar (within a given model) for high mass clusters 
and field galaxies, but significantly larger in low mass clusters. For 
comparison, the characteristic lumin osity of field gal axies in the 
K-band is Lk* = 6 x 10^°h~^LQ fcple et al.™200l'). The metal 
production in clusters and in the field is thus predicted to be domi- 
nated by galaxies with Lk ^ Lk* ■ 

In our models, the median galaxy luminosities and masses for 
metal injection into the ICM are almost identical to those for metal 
production. Our median stellar masses for metal production are 2-3 



times larger in the case of the most massive clusters than the median 
galaxy mass of I x lO^''ft.~^M0 f or metal (effectively oxygen) in- 
jection into the ICM calculated bv lDe Lucia. Kauffmarm & White I 
(2004J for a 10^^ h^^ AIq cluster. However, we are in agreement De 
Lucia et al. on the general conclusion that the ICM metal enrich- 
ment is dom inated by fair ly massive galaxies, in contrast to sugges- 
tions (e.g. Gar nett l2002h that metal ejection by dwarf galaxies has 
dominated the chemical enrichment of the ICM. 



8 SUMMARY 

We have investigated the metal enrichment of the ICM for clus- 
ters of different masses using a semi-analytical model of galaxy 
formation and evolution, GALFORM (Cole et al. 2000). We have 
implemented in this model a detailed treatment of metal produc- 
tion by both SNe la and SNe II, based on the predictions of stellar 
evolution models for both stellar lifetimes and yields of different 
heavy elements. This enables us to follow the evolution of differ- 
ent elements (O, Mg, Si and Fe), including the long time delays for 
SNe la. Our model consistently calculates the cycling of gas and 
metals between stars, cold gas and hot gas. In order to reproduce 
the observed present-day galaxy luminosity function using a realis- 
tic baryon fraction in the models, we have considered two alterna- 
tive mechanisms to prevent the overproduction of bright galaxies: 
(a) galactic superwinds which expel gas and metals from dark ha- 
los, with subsequent recapture in very massive halos, and (b) ther- 
mal conduction which transfers heat from the outer to the central 
regions of hot gas halos, thereby supressing gas cooling. We con- 
sidered two different IMFs for bursts of star formation: the Kenni- 
cutt IMF that we also adopt for quiescent disk star formation, and 
a top-heavy, flat IMF with x — 0. We have found that the form of 
the IMF is strongly constrained by the metal abundances of O, Mg, 
Si and Fe in the ICM. Our main conclusions can be simimarized as 
follows: 

(1) A top-heavy IMF in starbursts is required in order to obtain 
ICM metal abundances as high as those observed, when the IMF for 
quiescent star formation in disks is fixed to be a standard, bottom- 
heavy IMF, such as the Kennicutt IMF. A burst IMF with x = re- 
produces the observed metal abundances in the ICM. We have also 
tried the Arimoto-Yoshii IMF with a; = 1 in bursts, and confirmed 
that it fails to produce sufficient metals. Considering that the IMF 
for quiescent star formation has been w ell constrained by many 
studies, including a semi-analytical model I' Nagashima & Okamotd 
i2004) . this is strong evidence in support of a top-heavy IMF with 
a; ~ in starbursts. 

(2) Both the superwind and conduction models reproduce the 
observed ICM metal abundances and abundance ratios. The super- 
wind model predicts lower metal abundances in low temperature 
(T < 1 keV) clusters. Therefore if the drop in [Fe/H] for clusters 
with T < 1 keV found using A5CA is real, this could be evidence 
for the action of superwinds. The current data are still limited, how- 
ever, so future observations, particularly of low temperature clus- 
ters, would strongly constrain the superwind model. 

(3) Our model is in broad agreement with the observed ratio of 
ICM gas mass to the total 6j-band stellar luminosity in clusters, 
Afhot/ii),. In particular, the superwind model predicts a decrease 
of Mhot/ifcj in low temperature clusters, as is observed, because 
the superwind expels gas from these clusters. 

(4) The iron metallicity, [Fe/H], is predicted to decrease only 
slowly with increasing redshift. The difference in [Fe/H] between 
z — and 1 is at most 0.2 dex. This mild evolution is broadly 
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consistent with the recent observations by iTozzi et alJ i2003h . In 
contrast with [Fe/H], [0/H] is almost independent of redshift over 
this range, so [0/Fe] increases with redshift. 

(5) The metal production histories of clusters having the same 
present-day mass are very similar. In contrast, the cluster mass as- 
sembly histories show a wide variety. This suggests that the variety 
of metallicity gradients observed in clusters may be caused not only 
by observational errors but also by the variety of assembly histories. 
Presumably this is the origin of the diversity in the observed metal- 
licities of present-day clusters, because most observations measure 
metallicities only at cluster centres. 

(6) In clusters and in the field, metal production is dominated by 
galaxies with K-band luminosities similar to or brighter than L*. In 
low-mass clusters, the contribution of the central galaxy to metal 
production is very important, but this effect is weaker in high-mass 
clusters. 

Throughout this paper, we have neglected any contribution of 
intracluster stars to the metal enrichment of the ICM, since in our 
model all stars are in galaxies. Based on recent measurements of 
the intracluster light, it has been suggested that intracluster stars 
might have produced a large fraction of th e iron observed in the 
ICM tea ritskv. Gonzalez & Zabludofi2004) . The intracluster stars 
would presumably originate from material tidally stripped from 
galaxies or ejected during violent mergers, so the luminosity func- 
tion of galaxies in clusters would also be affected. In our model, 
tidal stripping of stars from galaxies would not affect the total mass 
of iron produced (if all other parameters remain the same), but it 
might affect the fraction of the iron which is ejected into the ICM. 
This possibility clearly deserves further observational and theoreti- 
cal study. 

It is also evident that, in order to better constrain our model, 
we need more accurate ICM metallicity measurements for larger 
samples of clusters. These measurements need to cover the outer 
parts of clusters as well as the central regions, in order derive ac- 
curate global metallicities. The observational situation will be sub- 
stantially improved by the Astro-E2 satellite scheduled for launch 
in 2005, as well as by further data from XMM-Newton and Chan- 
dra. Further constraints on the model can be obtained by examining 
the abundances of different heavy elements in galaxies, and we will 
investigate this latter topic in future papers. 
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APPENDIX A: CORRECTION OF THE OBSERVATIONAL 
DATA FOR METALLICITY GRADIENTS 

Some clusters are known to have substantial metallicity gradi- 
ents. The so-called 'cooling flow' clusters show enhancements 
in the metallicities of some elements toward the cluster cen- 
tre (e.g., Tamura et al. 2001; Matsushita, Finoguenov & Bo hrmgea 
^003; De Grandi et al. 2004). Tamura et al. 1 2001, 2004) have re- 
ported that in cooling flow clusters O and Mg among the a ele- 
ments do not show the metallicity gradients, whereas Fe and Si 
tend to show gradients, at least from the cluster centre to around 
half of the virial radius, despite the fact that Si is an a element. 
Non-cooling flow clusters generally do not show metallicity gra- 
dients for any elements. Our model predicts the global metallicity 
of the ICM, averaged over the whole cluster Therefore we need to 
correct the observed metallicities to obtain values averaged over the 
who le cluster, becaus e most observations give only central values. 

iDe Grandi et alj 12004.) compiled Fe metallicity gradients in 
12 cooling flow clusters. Based on their results, we assume that all 
cooling flow clusters have the same shape of metallicity gradients 
in both Fe and Si, 



Z{r) = Zo 



1+1^ 

Xc 



(Al) 



where x — r/r2oo, r2oo is the radius within which the mean over- 
density is 200 times the critical density, Xc = 0.04 and a — —0.18. 
We define the average metallicity within radius r as. 



Z{r)p,{rydr 



(A2) 



where pg (r) is the gas density. We assume that the observed metal- 
licities correspond to average values within some radius r, Z{^ r), 
and then use the assumed form of the metallicity profile (equation 
(eqn:metalgrad))to correct these to mean values Z{^ rvir) within 
the virial radius r-vir. The central metallicity Zq is thus treated as a 
parameter chosen so as to match an observed metallicity Z{^ r). 



Al De Grandi et al.'s sample 

iDe Grandi et alJ i2004b give the iron masses of individual clusters 
as a function of A, which is a spherical overdensity. 



3Mtot(< ta) 



(A3) 



where A/tot (< t) is the total mass within a radius r and 
Pc,z is the critical density of the universe at redshift z. \]s,- 
ing the data on the distributions o f total and gas mass given by 



lEttori. De Grandi & MolendiH2002h . we obtain an averaged metal- 
licity for each cluster within rA,Z{^ ta). The gas density dis- 
tribution, Pg (r) is assumed to follow an isothermal profile with a 
core. 
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To obtai n ta, we use the structural parameters given by 
pittori, De Grandi & Molendi (2002) assuming an NFW profile 
iNavarro. Frenk & Whit&.1997.) for the total density distribution. 
We have found that the use of structural parameters assuming a 
King profile hardly affects averaged metallicities. is assumed to 
be the same as the scale radius of the NFW profile, but the results 
are not sensitive to the exact value chosen. 

The metallicity at the cluster centre, Zq in eg. lAU . depends on 
the individual cluster. Using the above equations and the values at 
''2500, at which the observed data are given in the above papers, we 
obtain the average metallicity Z within the virial radius (r 5C rvir, 
corresponding to A ~ 100 in a ACDM universe), 

2 _ -^Fe,2500 Z{^ rvir) 
A^gas,2500 '■2500) 

The average [Fe/H] is easily estimated from Z. 



so the Fe-L line dominates in abundance estimates. Lower temper- 
(A4) ature gas emits more photons in the Fe-L line than in the Fe-K 
line. This has the consequence that for clusters with cooling cores, 
the observationally-estimated iron abundance is weighted towards 
the value in the low temperature component, that is, the central 
region, where the iron abundance is typically higher than the av- 
erage over the whole cluster. Thus, as Fukazawa et al. 1 1998) have 
done, the data for the central region should be removed to get a 
fairer estimate of the me tallicity averaged over the whole cluster. 
iBaumgartner et all i2003h have not done this, and this may be the 
reason why their data show an enhancement of the iron metallic- 
ity at intermediate cluster temperatures. For this reason, we do not 
include their data on Fe abundances in our comparison with the 
models. Their data on Si abundances seem to be less affected by 
this problem, so we do include those data in our comparison. 

(A5) jj^jg paper has been typeset from a TjiXy LTj3£ file prepared by the 
author. 



A2 Peterson et al.'s sample 

Peterson et al. (2003) observed only the central parts of clusters 
within a 1' angular radius. Their observed metallicities of Fe and 
Si are assumed to correspond to Z{^ r^/ ), where r^/ is the radius at 
which the angular size is equal to 1 arcmin. The average metallicity 
within the virial radius is then obtained from 

Z = ZobB-^— -• (AO) 

The parameters regarding redshift and structure required to obtain 
averaged metallicities are tabulated in their paper. 

Recent observations suggest that the metallicity gradients of 
Si are very similar to those of F e (Tamuraetal. 200 iL E004I : 
iMatsushita. Finoguenov & Bohringer.,200 3). Therefore we assume 
that Si has the same shape of metallicity profile as Fe. As already 
mentioned, because O and Mg seem not to have any metallicity 
gradients, we regard the observed values of central metallicities as 
equivalent to the metallicities averaged over the whole cluster. 



A3 Fukazawa et al.'s sample 

The measurements of Fe and Si abundances by Fukazawa et al. 
(1998) include data from the outer parts of the clusters. We there- 
fore assume that their abundance measurements are close to global 
averages, and do not apply any corrections. 



A4 Baumgartner et al.'s sample 

IBaumgartner et all i2003h obtained average abundances of several 
elements, including Fe and Si, in bins of cluster temperature by per- 
forming a stacking analysis of observations of 273 clusters from 
the A5CA archives. Although their cluster sample is large, metal- 
licity gradients combined with the large point spread function for 
A5CA cause problems for interpreting their results, especially for 
Fe. The iron abundance is usually estimated from the Fe-K or Fe- 
L lines. For high-temperature clusters, the photon flux in the Fe-K 
line is larger than that in the Fe-L line, and the Fe-K line yields 
accurate measurements of the iron abundance. On the other hand, 
for low and intermediate temperature clusters, the Fe-K line be- 
comes very weak, and the photon flux is larger in the Fe-L line. 



